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It is shown that nonadiabatic fluctuations of the sohton lat- 
tice in the spin-Peierls system CuGeOa lead to an important 
reduction of the NMR line widths. These fluctuations are the 
zero-point motion of the massless phasonic excitations. Fur- 
thermore, we show that the discrepancy of X-ray and NMR 
soliton widths can be understood as the difference between a 
distortive and a magnetic width. Their ratio is controlled by 
the frustration of the spin system. By this work, theoretical 
and experimental results can be reconciled in two important 
points. 

75.10.Jm, 75.30.Fv, 75.50.Ee, 75.30.Kz 



I. INTRODUCTION 

Already 20 years ago, the physics around the spin- 
Peierls transition fascinated many researchers (for a re- 
view, see Ref. [y). In particular the incomniensurably 
modulated I phase attracted much interest (e.g. [p[-p^). 
Yet detailed experimental investigations of the nature of 
this phase were not possible at that time. The first spin- 
Peierls transition in an inorganic compound, CuGeOa, 
was found only five years ago [H. This made a multi- 
tude of experimental investigations possible (for a review, 
see Ref. |5|). 

In particular, direct X-ray experiments in the I phase 
were performed by Kiryukhin and Keimer which permit- 
ted for the first time to detect the incommensurability of 
the distortion in fc-space |Q. Even more, it was possible 
to look at the structure of the soliton lattice modulation 
by measuring the intensity of the third harmonic ||l3,|l8[ . 

On the other side, Fagot-Revurat et al. |19] were able 
to measure the distribution of local magnetizations in 
CuGeOs in a beautiful NMR experiment. In a refined 
version it was now possible to deduce from such results 
the shape and the amplitude of the magnetic part of a 
soliton [ pOpl| ]. Three discrepancies to the conventional 
theories became apparent. The first concerns the ampli- 
tudes of the local magnetizations which is experimentally 
found to be much lower (factor 4 to 6) than predicted. 
Second, the X-ray soliton width (13.6 ± 0.3) is appre- 
ciably larger than the NMR soliton width ranging from 
6 to 10. Third, the widths are all larger than the ones 
theoretically predicted. 



In the present work we will solve the first two discrep- 
ancies and argue with Zang et al. |22| that the remaining 
problems are connected to the neglect of interchain cou- 
plings. 

To fix the diction let us state that we use the term 
soliton for the combination of a zero in the modulated 
distortion and the concomitant localized, bound spinon 
[0,|U23] . The distortive soliton width is the width of the 
kink-like zero of the modulated distortions. The mag- 
netic soliton width is the spatial width of the local mag- 
netizations |Q,||,|^,|3| . The incommensurate modulation 
in the I phase is viewed as an equidistant array (lattice) 
of solitons. 

The paper is set up as follows. In section II we dis- 
cuss the fact that the NMR line widths are much smaller 
than the theoretical ones and give an explanation for it 
in terms of the zero-point motion of the soliton lattice. 
Direct numerical calculations based on DMRG are shown 
in section III. They permit to address the second main 
point, namely the discrepancy between the soliton widths 
as measured by X-ray and by NMR. In section IV a de- 
tailed comparison to recent experimental data |^,|l| will 
be presented. The concluding section contains a discus- 
sion of the open questions, namely the role of interchain 
couplings, and a summary of our results. 



II. AVERAGING DUE TO PHASONS 

In Ref. |jl9] the results were interpreted by fitting 
them to a Hartree-Fock theory of Fujita and Machida 
]l2[ . Fujita and Machida did not take into account that 
the expectation values which occur in the Hartree-Fock 
self-consistency problem become non-uniform in a non- 
uniform phase [ p5[ . Thus they missed the important 
point that the local magnetizations rrii := (Sf) are so 
strongly alternating that they are even antiparallel on 
every second site to the applied external magnetic field. 
The amplitude of the alternating component is strongly 
enhanced compared to the XY model as was predicted in 
a number of investigations 0,E3-p^,E3l . 

Theoretically, there is no doubt that the spin- isotropic 
model has to be used to describe cuprate systems. Ex- 
perimentally, however, the amplitude of an XY model fits 
much better than the enhanced amplitude of the isotropic 
XYZ model. The discrepancy can be explained by the 



fluctuations of soliton lattice which are induced by the 
presence of the so-caUed phasons ||2^ . 

The importance of zero-point motion of the crystal lat- 
tice in Peierls systems was aheady noted by McKenzie 
and Wilkins |28]. In the present work we wiU focus on 
the zero-point motion of the phasons. The phasons arc 
very similar to phonons. If the deviation from commensu- 
rability of the wave vector characterizing the modulation 
d := Iq — 7r| is small (low soliton concentration) and the 
soliton width ^ is large a continuum approach can be used 
Pp^ . In this approach the discreteness of the underly- 
ing lattice does not appear. Thus the continuum Hamil- 
tonian is invariant under continuous translations along 
the chains. The incommensurate modulation breaks this 
continuous symmetry giving rise to Goldstone bosons, 
the so-called phasons p9]. They refer to oscillations of 
the solitons about the equilibrium positions in their lat- 
tice just as phonons refer to oscillations of the atoms 
about the equilibrium positions in the crystal lattice. 
There is one important difference between phonons and 
phasons. There exist in general three phonon branches 
(2 transversal, 1 longitudinal) corresponding to the three 
spatial dimensions into which an atom can be moved. 
But there is only one phason branch since the modula- 
tion can only be moved along the chains. Note that this 
does not concern the fact that the phasons have a non- 
degenerate dispersion uj(k) which depends on a three di- 
mensional vector k. Like phonons the dispersion is linear 
in k for small values of fc, i.e. 



^ ={cccK+Cyk'y + 2\c,\ki)/p 



(1) 



in the notation of Ref. p3]. Hence, the phasons give rise 
to a T^ contribution in the specific heat |^ which is 
indeed experimentally observed pffl . 

Before we proceed further we discuss briefly the effect 
of pinning. First we like to emphasize that there are two 
possible sources of pinning. The first one is pinning to 
the discrete lattice structure. The second is pinning to 
defects. 

The first mechanism enters since the continuous trans- 
lational invariance along the chains is given only in the 
continuum treatment which represents a certain approxi- 
mation. So the phasons are only quasi Goldstone bosons 
of a quasi continuous symmetry breaking. Yet treating 
the incommensurate modulations as continuously trans- 
lational invariant, i.e. shifting them without energy cost, 
is an excellent approximation if the solitonic width f is 
not too small. In the course of our previous calculations 
p4| we noted that the energy difference between a mod- 
ulation with the zero on a site and a modulation with the 
zero between two sites is of the order Jexp(— C^) where 
C is some constant of the order of the inverse lattice con- 
stant c~^ in chain direction. Hence, for ^ « 10c this 
energy difference becomes negligibly small. 

From the experimental point of view we come to the 
same conclusion. If there were a pinning of the modu- 
lation to the lattice structure the modulation would be 



commensurate with a period Lc where L is an integer. 
This would imply that at most L discrete local magne- 
tization values mi occur. The experimental resolution, 
however, is sufficiently high to exclude this scenario since 
not a number of isolated peaks but a continuous distri- 
bution is observed in the NMR response [£0|,^ . 

The second mechanism is pinning to defects which 
break the translational invariance. Such an effect is cer- 
tainly present but it is negligibly small in the pure sam- 
ples. To obtain an estimate we argue that the defect con- 
centration X corresponds to a typical distance between 
two defects of / = c/x where c is the lattice constant in 
chain direction. This means that phasons with a wave 
vector below k^ ~ 2tt/1 = 2ttx/c do not exist. This can 
be viewed as the effect of a gap Apin induced by defect 
pinning. From eq. (H) we obtain the estimate 
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where we used c? = uq^q^. Atrip = 'h\/unl P = 24K, 
i^z = 0.69nm H and c = 0.294nm |l|. An upper 
bound for the defect concentration x in the pure samples 
investigated is 10^'^. So the defect pinning gap Apjn is 
lower than 0.5K. We conclude that defect pinning will 
become important only below T k, 0.5K and does not 
need to be considered here. 

The soliton comprises a zero of the modulated distor- 
tions and a spinon bound to this zero [E3[. If there is a 
zero-point motion of the phasons this implies a certain 
motion of both the lattice distortions and the magnetic 
structure. It is plausible to assume that a certain aver- 
aging occurs which reduces the amplitude of the alter- 
nating magnetizations. This idea was first introduced in 
Ref. [p5| to explain the difference between observed and 
computed magnetization pattern. (Note that Kiryukhin 
et al. jlj] discussed a certain phasonic averaging linked 
to defects. Our approach does not rely on defects.) Here 
we present the detailed calculation and further estimates. 

Let us assume that the local magnetizations rrii can 
be described by two smoothly varying functions a{r) 
and u{r) which provide the alternating and the non- 
alternating component, respectively 



TTii = a(ri) cos(7rrj) -|- u{ri) 



(5) 



where we set the lattice constant to unity and r^ de- 
notes the component along the chains. The continuum 
approach results p^ , p2| are in fact of the form (0). Eq. 
(g) is the adiabatic result describing the completely static 
situation without phasons. Let us introduce now the 
phase variable 6(ri) where we use the hat to indicate 
that it is an operator as is the position of a harmonic 
oscillator. Thus (ph becomes 



rhi = a{ri) cos{nri + 9(ri)) + u{ri) 



(6) 



In principle, the shift Q{fi) has to be inserted in the 
functions a{ri) and u{ri) as well. But these functions 
are slowly varying so that the influence of the shift on 
them is negligible. Assuming furthermore that the NMR 
experiments measure on a relatively long time scale we 
conclude that the local magnetization to™^ seen in ex- 
periment is simply the expectation value 

m^'^P = (mi) (7) 

In the harmonic approximation the phase fluctuations 
are interactionless bosons and the operators &{fi) are 
linear combinations of the bosonic creation and annihi- 
lation operators. Then it is straightforward to compute 
the expectation value of the cosine 

(cos(7rrj -I- Q{n))) = cos(7rrj)(cos(e(r;))) (8a) 

= cos(7rr,) exp(-(e2(r,))/2) . (8b) 

Since the dominant fluctuations are those at long wave 
lengths the dependence of (O^(ri)) on the site index i 
should not be important. Hence we introduce 



exp 



^zmM)<^ 



(9) 



where TV denotes the number of sites in one chain. The 
reduction factor 7' is similar to a Debye- Waller factor 
which accounts for the non- vanishing atomic motion due 
to phonons. It reduces the amplitude of the alternating 
component only. From eqs. (PPB,0) we find 



"^T^ = I'a-in) cos(7rrj) -I- u(rj) 



(10) 



where the essential amendment compared to (0) is the re- 
duction of the alternating component by 7'. It is plausi- 
ble to explain the discrepancy between experimental and 
adiabatic theoretical amplitude by the zero-point motion 
of the phasons which leads to a finite value of (G)^(ri)) 
and hence to 7' < 1. We will present further support for 
this idea in section IV. 

Before turning to estimates for 7' we point out how 
one can take the reduction 7' into account if the result of 
the adiabatic calculation is not given in the form (|^) but 
as a set of discrete values {rrii}. This is the case for any 
direct adiabatic numerical treatment (see, e.g., section 
III) which does not use the continuum approach. Then 
one has to deduce in a first step estimates for the slowly 
varying functions a{ri) and u{ri) from the rrii. The most 
natural way to do this is by taking local averages 

a(ri) = mi/2 - (mi_i -|-mi+i)/4 (11a) 

-u(ri) = mi/2 -I- (mi_i -|-mi+i)/4 . (lib) 

After application of 7' to u{ri) as in ( |lO| ) one obtains 



:7'a(ri) +u{r^) 



(12a) 



= (1 + -f')m,/2 + (1 - 7')(m,_i + m,+i)/4 (12b) 
= (1 - 27)mj + 7(mi_i +mi+i) (12c) 



with 7 = (1 - 7')/4. Eq. ( |l2c| ) is at the basis of the 
averaging of adjacent sites which we used previously P5| . 
Forster et al. even average completely over two adjacent 
sites rrii -^ {'rrii + Wi+i)/2 P7| . 

Now we turn to the calculation of 7'. Expressing the 
expectation value in the exponent in momentum space 
yields 



3 k 

?i ^ 1 A 2 



(13a) 
(13b) 



where we used 9(fc) = J'h/(2Muj(k)){a^ -f at) with the 
mass M . The first term in the bracket in ( |l3h| ) stands for 
the zero-point motion since it survives even for T — > 0. 
The second term in the bracket is the bosonic occupation 
number at finite temperature. Rearranging the exponen- 
tials yields 



7' = exp(-i:)/2) 

hv r 1 

2M J Z{k) 
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2k^T / (27r)3 



(14a) 
(14b) 



where v stands for the volume per spin site. 

For a comparison with experimental data we can fo- 
cus on the low-temperature behaviour of (jlj). A close 
inspection of (^ reveals that D = Di + D2 + 0{T^) 
where Di is constant and D2 is of order T^. With the 
help of the input from Ref. Q Di and D2 are deter- 
mined in the appendix. One obtains 
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D. 
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trip 



2\/2 uq£,qv'^/3 
{k^Tf fT 



6V2^"3uoAtrip \T 



(15a) 
(15b) 



where Atrip is the singlet-triplet gap and ^0 ^^nd uq a 
characteristic length and characteristic energy per vol- 
ume, respectively, defined and given in Ref. [E9|. The 
value of Di for CuGeOs is 3.71 and the characteristic 
temperature T* is 16. 9K. 

The value of Z?i leads to a reduction of the alternating 
component 7' = 0.16 and the parameter 7 takes the value 
0.21. This is in very good agreement with the values 0.19 
and 0.20 which we found previously by fitting theory to 
experiment [ p5[ . 

Further support is gained from the estimate for T* . In 
Fig. n^ the temperature dependence of the widths W of 
the NMR lines are shown which are dominated by the 
alternating component. Hence they are expected to be 
proportional to 7' — exp(— (r/T*)^/2) to which they are 
compared. 




40.0 , 60.0 

FIG. 1. Solid line: theory from eq. ( |l5h| ); symbols: ex- 
perimental data (the error is less than a fifth of the symbol 
size). 

The small number of data points at low T does not 
allow for a complete quantitative test of the theoretical 
prediction. But the order of magnitude of the reduction 
of the line width at low temperatures is the correct one. 

If there were an adiabatic treatment of the incommen- 
surate phase at finite temperatures this would yield a line 
width which saturates exponentially at low temperatures 
due to the gap in the spin system p2| , p4| . Fig. [| indicates 
clearly that there is no exponential saturation of the line 
width for T ^ 0, but a behaviour W-W{T = 0) ex -T^. 
This is a direct evidence for the presence of low-lying, 
gapless fluctuations. The quadratic beh aviour in T of 
the decrease corroborates the result ( |l5b| ) and hence the 
existence of a three-dimensional nondegenerate disper- 
sion. 

At temperatures closer to the transition the line width 
is reduced much more strongly. This is due to the fluctu- 
ations in the spin system itself which are not within the 
scope of the present treatment. 

In spite of the crudeness of the present estimates the 
agreement gives evidence that the basic idea, phasonic 
fluctuations leading to an average of the local magneti- 
zations, is correct. 



III. NUMERICAL RESULTS 

Here we present some numerical results obtained by 
DMRG. The calculation treats the phonons adiabatically, 
i.e. the local distortions Si are real numbers which are 
found by minimizing the ground state energy of the fol- 
lowing Hamiltonian |26|j24|] 



^ = E 



J ((1 + (5i)SiSj+i + aSiSj 



K 



hSf 



(16) 

The Si are determined self-consistently [p6|,p3. Numer- 
ically, this is achieved by iteration. Since it is found 



that an array of equidistant solitons represents the ener- 
getically most favorable configuration |^,Q we use as 
initial distortion Si oc cos(gri) with q = tt + 2TTm, m 
being the average magnetization or Si oc sign(cos((7ri)). 
After about 10 iterations, one reaches a stable distortion 
pattern. This pattern shows no discernible differences 
depending on which start configuration has been used. 

For the calculation of the lowest state with S* = 1 for a 
given modulation we use the finite size algorithm p3,p5| . 
In each iteration we keep m = 64 states. Periodic bound- 
ary conditions are applied. Keeping m = 128 states leads 
to a change of the calculated energy of the order of 10~^. 
For large chain lengths and large dimerizations the en- 
ergy change due to a change of the distance by a few sites 
between two neighboring solitons is very small. This is 
due to the exponential localization of the solitons. So 
care has to be taken to avoid spurious shifts which hin- 
der the following fit analysis. 

In Fig. H the local distortions and the local magneti- 
zations are shown as they are found after the procedure 
described above. The solid lines are fits of the following 
form 
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(17a) 
(17b) 



where the parameters W , R, S, fcm, and fed are taken to 
be the fit parameters. The periodicity L of the solitons 
is fixed by the average magnetization m = 1/L which 
results in the relation defining the soliton widths ^m and 
^d as a function of L and fc,n/d 



4mfcj„/dii'(fcm/d)Cm/d 



1 



(18) 



Another relation can be deduced from the average of the 
dn- function which is t: /{2K{k)) [Ba] which in the contin- 
uum limit is related to the average magnetization. One 
obtains 



ttW 



mK{k^) 



(19) 



a relation which was well fulfilled (within 1 to 4%) by our 
fit parameters. 

The motivation for the equations (|l7| ) is twofold. First, 
such fits are used to describe the experimental data 
|19|-pi|. Second, the continuum approaches provide re- 
sults of the above kind [|J3[0,|,|l|,|2| . The continuum 
results predict that the magnetic and the distortive pa- 
rameters are identical 



Km — Kd <=> i;in — ^d 



(20) 



The DMRG results can be fitted very well by (jlTl). In 
this respect, the continuum approach works fine. But in 
other respects it fails. 
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FIG. 2. Upper panel: local distortions. Symbols stand for 
the self-consistent DMRG-result at K — 18J and a = 0.35; 



solid line steins from (17b) with S = 0.014, kd = 0.959, and 
^d = 10.5. Lower panel: local magnetizations; the solid line 
stems from ( p7a] ) with W = 0.21, R = 5.0, fcm = 0.992, and 
U = 7.9. 

Taking, for instance, the results for the amphtudes 
of the non-alternating and the alternating component 
T2l we get for W and R 



W 



R 



/ 2 Atrip 

Trws 
' 27rws 

Atrip 



(21a) 
(21b) 



Inserting some reasonable number for CuGeOa Atrip = 
24K, fc « 1, ws = f J(l - l-12a) H) with a = 0.35, J = 
160-?^ yields W — 0.32 and R = 6.3. So the agreement 
with the numerical results presented in Fig. [^ is not good, 
but the right order of magnitude is reproduced. Another 
indication that eqs. (pi] ) have to be extended is the fact 
that in a ID approach for a > Uc = 0.2142 a gap opens 
and vs is no longer well defined. 

Similar to the amplitudes, the unique soliton width is 
given in the framework of the present continuum theories 
by 



between about 10 and 6 with the higher number close to 
the transition. In fair accordance with our calculation, 
the X-ray measurement (susceptible to the distortion) 
yields a value about 1.4 larger than the value obtained 
by NMR [El| (susceptible to the local magnetizations). 
We will elucidate this issue further at the end of the fol- 
lowing section. 



IV. COMPARISON TO EXPERIMENTAL DATA 

In this section we attempt a quantitative comparison to 
the experimental results obtained recently in high qual- 
ity pl|. We will do this on the basis of the Hamilto- 
nian p6[), i.e. a magnetic one dimensional Hamiltonian 
with adiabatic phonon treatment. We use a = 0.35 
and K ~ 18J. The value of frustration results from 
fits of the temperature dependence of the susceptibility 
[^,Q. The value of K is then necessary to account 
for the amount of dimerization 6 ~ 0.014 which yields 
the correct size of the singlet-triplet gap in CuGeOa of 
Atrip/ J = 2'iK/160K = 0.15. This value of K provides 
also a reasonable estimate for the critical magnetic field 
at T = [^. One must be aware, however, that the in- 
clusion of higher dimensional magnetic couplings implies 
a considerably larger dimerization to keep the same gap 
pO| . A hint that the actual K value could be smaller 
is provided by the adiabatic analysis of the spin-Peierls 
temperature Tsp. There a value of A' « 11 had to be 
used to reproduce Tsp ~ 14.4 [|2|. We will discuss the 
effects of the neglect of the higher dimensional couplings 
in the concluding section. 

In Fig. H the amplit ude s computed numerically (filled 
circles) as defined in (17a) are contrasted to the exper- 
imental ones (filled squares; after p^). To obtain the 
computed values the self-consistently determined pat- 
ter ns a re least-square fitted with the functions given 
in (17a). The least square fits are very good (x^ ~ 
10^'^ — lO^"'), but not excellent. There are tiny devi- 
ations at the magnetic tails of the solitons, see Fi g. H . 
The fits for the lattice distortions on the basis of ( |l7b[ ) 
are better by two orders of magnitude. No deviations are 
discernible. 



S, = WS /Atrip 



(22) 



which takes the value 6.4 for the above numbers. Again, 
this is too low compared to our numerical results. 

The most important discrepancy is the difference be- 
tween the magnetic soliton width ^m = 7.9 and the dis- 
tortive soliton width ^d — 10.5. Note the ratio ^^d/Cm — 
1.33. This is very interesting because such a ratio can ex- 
plain the different experimental findings for the soliton 
width. By X-ray measurement [nTJ, ^d was determined 
to be 13.6 ± 0.3 whereas by NMR ^rn was found to vary 
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FIG. 3. Amplitudes as defined in eq. (17a). Filled circles: 
DMRG calculation for a = 0.35 and K = IS (corresponding 
to 5 = 0.014); open circles: by 7' = 0.25 renormalized DMRG 
data. Filled squares: experimental NMR data after Ref. | pi[ . 

For W and R, the striking difference between exper- 
iment and theory can be remedied by using a factor of 
7' = 0.25. The theoretical results renormalized in this 
way are shown as open circles in Fig. ^. The simulta- 
neous reduction of W and of R leaves the homogeneous 
part u{ri) of the local magnetization rrii unchanged but 
reduces the alternating part a{ri). This is exactly what 
we proposed in section II to be the effect of phasonic 
zero-point motion. The agreement between renormal- 
ized theory and experiment is good and provides further 
support for the phasonic averaging. 

In Fig. the results for the distortive and the mag- 
netic soliton widths are plotted. It is obvious that they 
are not the same as was assumed hitherto. Their ratio is 
to very good approximation constant for small magneti- 
zations but grows for larger m. This aspect is in quali- 
tative agreement with the experimental situation where 
there is also a striking difference between the X-ray result 
(distortion) and the NMR result (local magnetizations). 
The fact that the agreement is quantitatively not better 
can be attributed to several circumstances. The analysis 
of the X-ray data [[l7| was done on the assumption of con- 
stant soliton widths which is not justified. Furthermore, 
our theoretical analysis is still based on a one-dimensional 
model only. It is also obvious that the qualitative evolu- 
tion of the magnetic soliton width is not yet understood. 
Experimentally, it decreases with increasing magnetiza- 
tion whereas it increases in our computation. 
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FIG. 4. Soliton widths: ^d distortive; ^m magnetic, see eq. 
([l7|). Circles: DMRG results; Open squares: experimental 
NMR data after Ref. |2l|, filled square: value proposed in 
Ref. pj]. Upper panel: filled symbols are based on the dis- 
tortion; open symbols are based on the local magnetizations 
nii. Lower panel: ratios. 

Based on the lower panel of Fig. ^ we conclude that 
the ratio ^d/'Cm becomes constant if the solitons are suf- 
ficiently separated. Hence the ratio does not depend on 
TO or ii" as long as Xjm is large enough compared to ^d- 
Recall that solitons are exponentially localized objects. 
The amplitudes and spatial shapes do not depend on the 
overall energy scale J either. Thus the only parameter 
left is the frustration a. So we are led to an analysis of 
the ratio ^d/Cm as function of a. The results are depicted 
in Fig. I 
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FIG. 5. Dependence of the ratio Cd/Cm on the frustra- 
tion a. The values K are chosen such that the ^d ranges 
between 5 and 10, i.e. not too small but always way be- 
low the sample length oi L — 108. The values (a, K) are 
(0, 3), (0.1, 4), (0.2, 6.2), (0.241, 8), (0.35, 18), (0.5, 37). Solid 
line: regression with intercept 1.06 and slope 0.77. 



It is evident that the frustration is the important con- 
trol parameter for the ratio of ^d and ^m- The depen- 
dence seems to be roughly linear. It should be noted that 
even at zero frustration the two soliton widths are not 
identical. The experimental finding of a relatively large 
values of ^d/Cm is a strong indication for an important 
frustration in agreement with the x(T) fits |3^,Q. A 
quantitative analysis appears presently premature since 
the effects of higher dimensional couplings are not known 
yet. 

The deviation in the amplitudes as given in eqs. (EL) 
as well as the non-equality of the soliton widths suggests 
an extension of the existing continuum theories. The fact 
that the non-equality is strongly enhanced by the frustra- 
tion points clearly into the direction that the backscatter- 
ing cosine terms in the sine-Gordon models P;p2| is the 
origin of this effect. Analysing the self-consistent har- 
monic treatment of Nakano and Fukuyama Hg reveals 
one aspect where a difference between the magnetic and 
the distortive soliton can appear. This is the so far ne- 
glected spatial dependence of the renormalizing factors 
exp(-(e^)/2) and exp{-2{e'^)), respectively. (Note that 
this operator Q is not related to phasons of distortive 
origin. But it stands for gaussian magnetic fluctuations 
and their formal treatment bears many analogies to the 
treatment of phasons.) Further work along these lines is 
in progress. 



V. DISCUSSION 

In this work we calculated the effect of phasonic zero- 
point motion on the local magnetizations rrii. This nona- 
diabatic effect can be viewed as virtual oscillations about 
the static situation. It leads to a reduction of the alter- 
nating component of the rrii . The reduction factor 7' as 
in eq. (|^) is similar to a Debye- Waller factor which is 
induced by phononic motion. 

By the inclusion of the nonadiabatic phason motion 
we could explain the so far not understood difference be- 
tween experimental NMR line sha pes E9H 21t| and previ- 
ous theoretical predictions ||13| , p5| " P7y23[ . Hence the re- 
duction of the NMR line width compared to the adiabat- 
ically predicted one is a strong evidence for the impor- 
tance of nonadiabaticity in CuGeOs. Nonadiabaticity is 
also considered to be important for temperature depen- 
dent couplings and induced frustration (cf. 41-|43|). 

The second main point of our work is the difference 
between the distortive and the magnetic soliton width. 
To our knowledge, such a difference has so far not been 
reported in the literature. But it helps to understand 
the differing X-ray and NMR results. The ratio ^d/Cm 
is controlled by the frustration a. It is fairly close but 
not equal to unity at zero frustration and grows roughly 
linear to almost 1.5 at the Majumdar- Ghosh point a — 
1/2. Thus, further evidence for a sizable frustration in 
CuGeOa is provided. 



Aiming for quantitative agreement with experiment 
some questions remain. Choosing a small dimerization 
(the corresponding large K, respectively), which repro- 
duces the correct gap value and critical magnetic field in 
a one-dimensional model, leads to amplitudes W and R 
which have to be renormalized by about 7' = 0.25 and 
not by 7' — 0.16 as estimated in section II. Thus these 
amplitudes seem to be too small. The magnetic soliton 
width ^m, is too low by about 20%. 

Choosing, however, a larger dimerization (the corre- 
sponding small K, respectively), which reproduces the 
gap value resulting from averaging the dispersion u}{k) 
perpendicular to the chains along fcb [Q , leads to larger 
amplitudes W and R. These larger amplitudes can be 
nicely reconciled with the experimental findings by a 
value very close to 7' = 0.16. This is essentially what 
was done in the previous work p5[ . The argument in fa- 
vor of the larger gap is that the neglect of the interchain 
coupling would lead to this larger gap [^ . But the cor- 
responding soliton widths would be even lower (see eq. 
(p2|)) and hence in worse agreement with experiment. We 
conclude from this that it is not just a question of fine- 
tuning the constants to achieve agreement. The purely 
d = 1 model is not sufficient. 

The magnetic interchain couplings have to be taken 
into account as has become already apparent from the 
two ways to choose the relevant gap. The magnetic in- 
terchain coupling has the two following effects. First, 
it reduces the observed gap value for a given dimeriza- 
tion |2^ . Thus a larger dimerization (smaller K) should 
be used. This is consistent with a number of recent in- 
vestigations Q,Q. The larger gap will lead to larger 
amplitudes. Second, Zang et al. |22| argued that the in- 
terchain coupling enhances the soliton widths which is 
exactly what is needed to reconcile theory and experi- 
ment. This enhancement is qualitatively easy to under- 
stand since an interchain coupling will favor the appear- 
ance of staggered magnetization. Regions of staggered 
magnetization are those were the solitons are, see Fig. 
0. So these regions are extended which implies an en- 
hancement of the soliton width. Hence, we come to the 
conclusion that a larger dimerization (smaller K) plus an 
appropriate treatment of the interchain couplings should 
yield results in quantitative agreement with experiment. 
There are also preliminary results indicating that the in- 
clusion of interchain coupling changes the increase of ^m 
as function of m to a decrease which is what is experi- 
mentally seen (cf. Fig. ||). 

Summarizing, two main points are shown in the 
present work. The first is the importance of the pha- 
sonic fluctuations providing a natural explanation for the 
large amplitude differences between theory and experi- 
ment. The second is the difference between distortive 
(^d) and magnetic (Cm) soliton width leading to a ratio 
^d/'Cm between 1.05 and 1.45 depending on frustration. 
Further work to include interchain coupling is called for. 
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APPENDIX: 



Based on eq. (14b) and Ref. ||29[ we present here es- 
timates for the two leading contributions Di and I?2- 
Focussing on the low-lying excitations we adopt (0) even 
though we do not know whether at larger values of k the 
dispersion is still described by (n^). 

The mass density p is proportional to Mjv. If each spin 
site moved in the same way on shifting the modulation 
M — vp were reasonable. Since, however, the modulation 
does not have the same amplitude for all sites an effective 
reduction has to be taken into account. For sinusoidal 
modulation this factor is 1/2, i.e. M = vp/2, due to 
the average value of sin^ . For all other modulations the 
reduction factor will be between 1/2 and 1; we choose 
1/2 since the higher harmonic content of the modulation 
is very small [|l7| . 

With k,j: -^ kx^/p/coo, ky -^ ky^/p/cy, and k^ -^ 
k^^yp/\2cz\ as well as a = roTo^g^ and fo = 
{(.Ox^Oy^OzY^^ we obtain 



D^^- 



o3/2 



PV^' 



■x^y^z 



d^k 1 



pV2 VCo^oTo 



(27r 
3/2 



3-COth(?ifc/(2fcBT)) 



(fik 1 



(27r)3 k 



coth 



fik \ 
2kBT) 



(Ala) 



(Alb) 



For the value of uq = ^'o^o and other values see Ref. [E9J. 
We split D in the zero temperature contribution Di and 
the temperature dependent rest D2 by using coth(a;/2) = 
l + 2/(exp(a;)-l). 



We estimate Di, the zero temperature reduction due 
to the zero-point motion of the phasons. It is difficult to 
compute Di rehably since, in principle, information of all 
the phasons is required, not only the lowest lying ones. 
But an estimate which determines the order of magnitude 
is possible by using an upper cutoff fcmax as in the Debye 
model of phonons M] . We start with 



Di^-J 77TZTI— -7 (A2a) 



h r d^k 1 

7) J WfZ{k) 

272^2 (^^^Y Jo 

hy/p 



r-fcmaj 



4V2^2 (^-^ 



/ kdk (A2b) 

a' "-max • (AzCJ 



The value for fcmax after rescaling is obtained from the 
number of possible phasonic excitations 

1 = w / j—:^ (A3a) 

(A3b) 




(A3c) 



Eqs. (A2c 



A3c) and the singlet-triplet gap Atrip 



hy^uo/p P9| together lead finally to (15a). It should 



be noted that the volume per spin site can be found from 
the lattice constants by w = abc/2 with a = 4.79- 10^^°7tt-, 
b = 8.40 • 10-i°m, and c = 2.94 • 10-^°™ [|l|. The divi- 
sion by 2 is necessary since there are two Cu ions in each 
unit cell. 

Next we compute Di where we are only interested in 
the leading temperature dependence. Thus the result 
depends only on the lowest lying excitations and the as- 
sumption of a dispersion as (|l|) is well justified. We ffiid 

, 3/2 ^oo/ ^^ 



^-71^ life j U-fe)-^j^^MA4a) 



= ^m^[w-) r(exp(fc)-l)-ifcdfc (A4b) 
sn-K^hp K^uqJ Jo 

_ VPJkBTf ^ jkBTf 

6V2n^lul^^ 6V2^3uoAtrip 



T' 



(A4c) 
(A4d) 



With the numbers Atrip = 24K, wq = 650mJ/cm^ and 
^0 = 0.31nm the characteristic temperature T* can be 
estimated to be 16.9 K. 



